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vacuum polarization function, Imn(t). This allows us to improve the corre- 
sponding contributions to the muon (or any other lepton) g — 2 anomaly and 
to the running QED constant on the Z, (Xqed(Mz)- This decreases the error 
in the evaluations by a factor between two and six for the high energy con- 
tribution, and by some 50% for the overall result. We find for the hadronic 

contributions a h = 6993.4 ± 110.0 x 10" 11 and Aa h = 272.59 ± 4.09 x 10" 4 . 
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I. INTRODUCTION 

There has been recently a renewed interest on two basic quantities, both related to 
the photon vacuum polarization: the hadronic contribution to the anomalous moment of the 
muon, a M , and the running QED charge a defined at the Z particle pole. In what respects 
the first, the reason is the advanced status of the experiment planned at Brookhaven which 
should improve the present accuracy by more than one order of magnitude; in what respects 
the second, because the quantity aQED(M^) plays a leading role in precision determination 
of electroweak parameters. 

In this paper we consider the contributions of the hadronic part of the photon vacuum 
polarization tensor to and a7Q££)(M§), that we wr he respectively as 

a( 2 )(had) , AagL(had , Mf) . (1.1) 

These quantities may be expressed in terms of the photon hadron vacuum polarization 

function 11^. One can then write a dispersion relation for this function in such a way that 

the experimentally accessible quantity 

= - 7* ^ hadrons) 

a(e + e~ — > /i + /U _ ) 

appears instead of 11^: 

roo 

a h = ajf^had) = / dt R(t) K(t) , (1.3) 
- "®°0* , Mi) = £ * A • (") 

Here, 

Kit) = ^ Cdz - Z } l ~ Z \ - . (1.5) 
V ; 3n 2 t Jo z 2 + (l- z)t/m* V ; 

The reason why we think we can improve on existing estimates of ah, Ao^ is that 
we may use the very reliable theoretical QCD calculations^ that have been extended with 



x The fact that QCD calculations are more precise than the use of experimental data was already 
noted in e.g. Ref. f2j. 

2 



great precision to regions previously unaccessible; in particular o to the regions right above 
threshold for heavy quark production (cc, bb) as well as the very low (t 1 / 2 ~ 1.2 GeV) 
energies. The extensions have been possible because of the high orders attained by the 
QCD calculations H and the use of the precisely known value of a s on the r mass H, 

a s (m 2 T ) = 0.33 ±0.05 . (1.6) 

Likewise, the use of the value ® of a s at the Z, 

a„{M%) = 0.119 ±0.003 , (1.7) 

allows us to reduce still further the uncertainties in the high energy regions. 

II. CONTRIBUTION FROM THE REGION t 1 / 2 < 1.1 GeV 

In the region t 1 / 2 < 1.1 GeV, perturbative QCD is clearly invalid. One thus has to rely 
on experiment, supplemented by old-fashioned hadron theory. The region may be further 
split into the p resonance region (say, t 1 / 2 < 0.8 GeV) and the rest. For we have, 

a h (t 1/2 < 0.8) = (4821 ± 24 ± 27) x 10~ n + (25 ± 3) x 10~ n . (2.1) 

Here the second term in the r.h.s. is the uj — p interference contribution. The first error 
is statistical; the second (when given) will be the systematic one. Eq. (2.1) presents the 
value reported in Ref. ||; other authors give compatible estimates. We will improve on (2.1) 
slightly later on. 

To this one has to add the contribution of the region 0.8 < t 1 / 2 < 1.1 which gives, using 
experimental data only, 

^(0.8 < t 1/2 < 1.1) = (1100 ± 94) x 10~ n . (2.2) 

This region, about which little can be done at present, presents the largest source of error. 
Combining Eqs. (2.1) and (2.2) we have, 
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(t 1/2 < 1.1 GeV) = (5947 ± 97 ± 27) x 1CT 11 . 



(2.3) 



For Acth, we have, taking the analysis of the first article in Ref. M 



Aa h (t 1/2 < 1.1 GeV) = (35 ± 0.3 ± 0.7) x 10" 



(2.4) 



The contribution of the region 1.1 < t 1 ! 2 < 2 GeV is the topic of next section. 

III. CONTRIBUTION FROM THE REGION 1.1 < t 1 / 2 < 2 GeV 

From a theoretical point of view, the contribution of this region presents a challenge, 
because one cannot use perturbative QCD reliably, and also because old-fashioned hadron 
theory does not describe very well the average experimental data points. In this section, 
we explain how one can combine these two theories and give reliable theoretical estimates, 
which can then be compared with previous estimates based on experimental data only : 



< xp (l.l < t 1 ' 2 < V2 GeV) = (278 ± 25) x lO" 11 , 
Aa^ xp (l.l < t 1/2 < y/2 GeV) = (13 ± 0.15 ± 0.8) x 10" 



(3.1) 
(3.2) 



Let us first estimate the contribution of this region to R(t) using perturbative QCD only: 
R(t) may be split as a sum over light quark flavors: 



R(t)= £ R q (t) , 

q=u,d,s 

and one may use the very precise high energy result for R q (t): 



(3.3) 



R^-(t) _= N c Q 2 ll-^l-v) 2 + ^l-v) :i + 



2 V 1 - v 



--)(l-v) 



3 9/-, 

- + -(1 - v) 

4 2 V ; 

C„a s 



C F a s 



71 



+ r 2 — + r 3 — 

7T V 7T / \ 7T 



,7 ('^V + 45.3 

IX J \ 7T 



(3.4) 



Here, 



4 



v = ^l-4m 2 (t)/t , ra = 1.986- 0.115 n/ , C F = 4/3 , 

f 3 = -6.637 -1.2^-0.005^- 1.24 (j^Qf] , (3.5) 

and m(t) and a s = a s (t) are the running mass and coupling constant which we may take 
to two loops i. The justification of this approximation lies in the fact that the light quarks 
are relativistic at these energies. Numerically, we take m u = m d = 0, and use @ 

m s (l GeV 2 ) = 0.19 GeV . (3.6) 

The result Rqcd(P) obtained in this approximation is plotted in Fig. 1. together with some 
experimental points. Although the experimental errors are large, it is clear that our curve 
Rqcd will not give a reliable estimate. Before we explain how this can be improved, let us 
first describe another approximation for R(t) based on old-fashioned hadron theory. 
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Fig. 1. 

To do this, we evaluate the contribution of the individual channels, in the quasi-two- 
body approximation. Thus, we have the two-body channels 7r + 7T K + K~, K°K , and the 
quasi-two-body channels p + 7r~, p~7r+, p°7r°, p°rj, ujti , <f)°TC°, (p°r), K*°K°, K* + K~, K*-R + . 

5 



The contribution of the first is estimated using the corresponding form factors, p or 
dominated: we obtain 



AjR(t) rcsonances = 1 L _ 4m|\ + | F ^|2 1 A 4m|A ? ( g ^ 



where 



m 2 p + m p m^ T p 



rri- \ — t — im p (t/ '4 — m 2 ) 3 / 2 £ ' (m 2 /4 — m 2 ) 3 / 2 



P m 2 + m^rriK^ T 4 

*K\t)- —2 — - — ——rn~A — ZTvurr > ?i = 



m 2 - t - i m^t/A - m 2 ,) 3 / 2 & ' Si (mJ/4 - m 2 ,) 3 / 2 ' 

Note that we have combined the contribution coming from both K + K~ and K°K° into the 
function Fk ■ 

The quasi- two-body channels contribution to R(t) is estimated in the vector-meson dom- 
inance approximation: 

AR(tr—^= £ Jf } . (3.8) 

Using the phenomenological interactions 

>~vp = 9vp e^F^V^p , (3.9) 

one can relate this to the radiative decay widths, 

U ^-,...«^K-4) 3 tV 2 (t_4m 2 ) 3 / 2 -^ Uj 

The decay widths are taken from experiment. The result -R(i)two-body we obtain in this 
approximation is also plotted in Fig. 1. It unfortunately does not represent an average of the 
experimental points. To bring these approximations closer to experiment, we construct the 
following interpolation curve which reproduces both evaluations in their ranges of validity: 
we note note that QCD works at the higher energies, and the old-fashioned hadron evaluation 
at the lower energies. Then, 

Rto\t) = ^two-body exp(-t 5 /V5t 5 ) + R QCD (t) [l-exp(-t 5 /v / 5^)] , 

t\ /2 = 1.25 , 1.3 , 1.35 GeV . (3.11) 
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The point to is chosen to fit experiment. These curves are drawn in Fig. 1 (solid curves) and 

represent a significant improvement from what we had previously. The numerical results 
obtained with these interpolation functions are: 

a h (l.l < t 1/2 < V2 GeV) = (265 ± 22) x KT 11 , (3.12) 

Aa h (l.l < t 1/2 < V2 GeV) = (15.3 ± 1.4) x 10" 4 . (3.13) 

1/2 

The central values have been obtained using R^o =1 -3GeV)^ anc j systematic errors have 



been estimated from the difference 



. (These errors are in 



^>(tj /2 =1.3 GeV) _ ^(tj /2 =1.35 GeV) 

our opinion over-estimated and should in principle be reduced. We will however not do that 
here). Comparing the numbers given in Eqs. (3.12, 3.13) with those given in Eqs. (3.1, 3.2) 
we see they agree, within errors. We however believe our evaluation to be more reliable than 
those estimates which were obtained using experimental data only, because of the known 
existence of large systematic errors of the last. 



IV. CONTRIBUTION FROM THE REGION t > 2 GeV 2 

For t > 2 GeV 2 , we distinguish two different contributions. The first is due to the J/ip 
and T bound states and the results can be taken from the first article of Ref. 0: 



a h ( J /if,) = (86 ± 4.1 ± 4) x KT 11 , (4.1) 

a h {T) = (1±0±0.1) x KT 11 , (4.2) 

Aa h (J/4>) = (11.34 ± 0.55 ± 0.61) x 10" 4 , (4.3) 

Aa h (T) = (1.18 ± 0.05 ± 0.06) x 10" 4 . (4.4) 

The second contribution comes from the continuum regions, i.e. the regions above qq 
thresholds. As before, we split R(t) as 

R(t)= J2 R i(t) ■ ( 4 - 5 ) 

q=u,d,s,c,b,t 
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For the light quarks u, d, s, only the high energy region needs to be considered; the 
corresponding expression for R q was given in Eq. (3.4). For the heavy quarks, c, b, t, we 
also need the value of R q for small and intermediate velocities. At low energy, the value of 
R q is known with great precision ®: 

2n{a s G 2 )\ 



+ 



v . J v(3-V*) , ( GV 37TV 2 ' 

IT 4 

TiC F a s /v 



E?{t) = N C Q 



C F a s }[l 



192 m 4 v e 



x[l + 2c (t)] 



1 - Am 2 /t 

93 - 10 n/ 
36 



I _ g-nCpOs/v ' 

a s = a s (t) 



(4.6) 



ai ±%A/2 ' 
1 H a s (£J 



7T 



A = 11 - 3™/ 



where 



c (*) 



In—- -1-2 7b 



(/ca) 2 (fca) 
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40 



+ 



= mt) , a = 2/mC F a s (t) . 



An exact expression for co(t) (which represents the radiative corrections) may be found 
in Ref. ||. Eq. (4.6) includes also nonperturbative corrections with (a s G 2 ) = 0.042 ± 
0.020 GeV 4 , and the evaluation is valid until these are of order unity, i.e., down to a critical 
velocity 



'2tt(« s G 2 )\ 1/6 



(4.7) 



v 192 4 m 4 / 

For the heavy quarks, we find it more convenient to reexpress the m in terms of the pole 
masses m: one has % 



a s (t) 




a s (m 2 )_ 





m(t) = m 



„ _ A 7o - A7i 
A = , 7o 



l Aa s {t) - (C F -A)a s (m 2 ) 



(4.8) 



4' 



-4, 4 = 102-yn/ 



7i 



202 20 
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with the values of the pole masses given by 



m c = 1570 ± 60 MeV , m b = 4906 ± 85 MeV 
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(4.9) 



and we take @ 

m t = 180 ± 12 GeV . (4.10) 

A crude estimate can now be obtained if one uses for R q (t) (q = c,b,t) the following ap- 
proximation: 

m/N \ R x t if v cnt <v<l/2 
R<p{t) = { q (4.11) 

I i?J e - if v > 1/2 

The expression is however discontinuous at v — 1/2. To obtain a smooth joining of the low 
and high energy regions, we can use the following two interpolation curves for R q (t): 

R q n \t) = {l-v)R^+vR h q e - , (4.12) 
R q IH \t) = (1 - v^R^- + v 3 R h q e - ■ (4.13) 

We consider Rg 11 ^ to be a very reasonable approximation, and take the difference between 
Ri 1 ' 11 ' 111 ) as an estimate of the systematic error due interpolations. We also take into account 
the error due to the input parameters a s (m T ), a s (Mz), and m t , as well as the error due to 
the non-perturbative contribution. The last is obtained by setting {a s G 2 ) = in Eq. (4.6) 
and allowing the velocity v of the heavy quarks to go to zero. Our best estimate for the 
contribution of the regions in the continuum is: 

^(continuum, t 1/2 > V2 GeV) = (826.4 ± 7.8 ± 9.8 ± 8) x 10~ n , (4.14) 
Aa h (continuum, t 1/2 > V2 GeV) = (226.6 ± 3.7 ± 1 ± 1) x 10" 4 . (4.15) 

The first systematic error is due to interpolations, the second is due to the error in the input 
parameters, and the last is the non-perturbative one. 

V. CONCLUSION 



Besides using our theoretical estimates for the high energy regions 

t 1/2 > 1.1 GeV, v^GeV 
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(5.1) 



for Imil(t) an extra (slight) improvement may be incorporated if we repeat the dispersive 
analysis of Ref. || in the p region using recent data. The reason for this improvement is the 
following: in that paper the p contribution to (say) a/ l (t 1 ^ 2 < 0.9 GeV), had been estimated 
from fits to the pion form factor F n alone (but both in the spacelike and timelike regions) 
or imposing also the values of the tttt phase shifts. The first method gave m p = 768 MeV, 
the second m p = 778 ± 2 MeV, both compatible (at the 2a level) with the then preferred 
experimental value m p = 769 ± 3 MeV. Although, as explained in Ref. f2|, the first method 
was considered more reliable, both were combined taking the difference between the two 
determinations as a measure of the systematic error of the calculation. Since the presently 
accepted experimental figure, m p = 768.1 ± 0.5 MeV clearly discriminates in favor of the 
method based on F n only, we can dispense with the (poorly known) tttc phase shifts and 
avoid the systematic error. We need only alter the evaluation by taking into account the 
change in the accepted value for the p width, T p = 151.5 ± 1.2 MeV from that used in Ref. 
0, 158 MeV (compatible with the 1985 experimental value, T p = 154 ±5 MeV but not with 
the presently preferred one). To first order this is easily taken into account as a variation 
AT/nm p = -0.27%. 

The results are summarized in the following tables where we also report, for purposes 
of comparison, the evaluations of the first (EJ) and third (MZ) papers of Ref. |l|]. We have 
composed quadratically statistical and systematic errors: there are so many of the last, and 
of such varied origins, that they may be taken to behave statistically on the average. For 
the muon anomaly we have, 

a h (t 1/2 <V2) = 6200.0 ± 101.4 x 10" 11 (This work) 
a h (t 1/2 <V2) = 6342.3 ± 137.7 x 10 -11 (EJ, Ref. [1]) 

and 

a h {t l/2 > V2) = 913.4 ± 14.9 x 10~ n (This work) 
a h (t 1/2 > V2) = 908.2 ± 77.1 x 10~ n (EJ, Ref. [1]) , 
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with the overall results 

a h = 7113.4 ± 102.5 x KT 11 (This work) 
a h = 7250.4 ± 157.6 x 10~ n (EJ, Ref. [1]) . 

which agree within the quoted errors. 

For the running QED coupling, one has 

Aa h (t 1/2 <V2) = 50.3 ± 1.6 x 10~ 4 (This work) 

Aa h (t 1/2 <V2)= 47.92 ± 1.06 x 10~ 4 (EJ, Ref. [1]) , 

and 

Aa h {t l/2 > V2) = 239.12 ± 4.1 x 10~ 4 (This work) 

Aa h (t 1/2 > V2) = 232.45 ± 6.45 x 10~ 4 (EJ, Ref. [1]) , 

and now the full results are 

Aa h = 289.42 ± 4.35 x 10~ 4 (This work) 
Aa h = 280.37 ± 6.54 x 10~ 4 (EJ, Ref. [1]) . 

Our central value for Aah is slightly higher than that of Ref. JI| (EJ), but deviates at the 
3-4 a— level from that given by the third article of Ref. |Ij]: 

Aa h = 273.2 ± 4.2 x 10" 4 (MZ, Ref. [1]) . 

We conclude this paper with a few comments on possible ways to improve the results. 
Certainly, better knowledge of a s both on the r and Z masses would increase the precision 
of the high energy evaluations; but most of the error comes from the region 0.8 < t 1 / 2 < s/2 
GeV. Even a modest improvement of a factor two in the error in the region between 0.8 
and 1.1 GeV would result in a substantial decrease, roughly by the same amount, of the 
overall error for a^, and about 20% for A«qed- This emphasizes the interest of some of 
the accelerators, projected or in construction, with the capability to explore these energy 
ranges. 
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